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Within the context of non-extensive thermostatistics, we use the factorization approxi- 
mation to study a recently proposed early universe test. A very restrictive bound upon the 
non-extensive parameter is presented: \q — 1| < 4.01 x 10~ 3 . 
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Although nonextensive formalisms are much in vogue in Physics, a very interesting one, the so-called Tsallis 
Thermostatistics (TT), seems to be one of the most actively studied. The proposed formalism basically relies 
upon two postulates 

• A new definition of a nonextensive entropy, 

w 

S„ = fc(l-X>?)/(flf-l). 

i 

• A new definition of expectation value, 

i 

TT introduces a new parameter, q £ 3?, which is usually called the nonextensivity parameter or the Tsallis 
q-index, and it contains the standard, extensive Boltzmann-Gibbs statistics, as a special case where q is 
taken to be unity. Up to present days, TT has been found to admit generalizations of some of the important 
concepts of statistical physics || , and to yield results which can explain some observational and experimental 
data where standard statistics is known to fail Q. All these efforts accelerate new attempts to find the 
physical meaning of the nonextensivity parameter, a long standing puzzle which has only now started to be 
clarified. On one side, some works have been devoted to the study of dynamical and dissipative systems 
On the other, an entirely new fractal canonical ensemble was introduced in order to relate TT with a scale 
invariant thermodynamics j(| . Of course, an alternative way to search for the meaning of q is related with the 
estimation of bounds in measurable physical systems. We mention the study of the microwave background 
radiation |7]|| , the Stefan-Boltzmann constant |],[l(J , the early Universe pj] , p^ | and the primordial neutron 
to baryon ratio in a cosmological expanding background In these applications, except in Refs. [8|JlO|], 
the quantal distribution functions of TT were obtained by an asymptotic approach of the kind — q) — > 0, 
where (3 is the usual inverse temperature. However, quantal distribution functions have been previously 
generalized by Buyukkihc et al. |14| using a rather straightforward procedure, referred to as Factorization 
Approximation. The formulae so obtained have proved to be simpler and more general than the ones 
derived with the common Tsallis et al. method [Q. Simpler, because the steps followed in the derivation 
are completely the same when compared with standard textbooks. More general, because the factorization 
5^ ' approach does not need a value of q ~ 1. 

Quantal distribution functions within the factorization approximation were, however, regarded as a rather 
rough technique because of a work by Pennini et al. [|l5| . They considered fermion and boson systems with 
very small occupation numbers. However, very recently, Wang and Le Mehaute analysed the problem 
in detail and showed that there exist a temperature interval, a forbidden zone, where the deviation from the 
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exact result is significant, but that outside this zone, the factorization approach results could be used with 
confidence. In addition, they verified that the magnitude of the forbidden zone remained constant with the 
increase of the number of particles, contrary to the result stated in Ref. jl5). This fact might motivate new 
efforts for the study of macroscopic systems (where the number of particles is ~ 10 23 ) within the simpler 
approach. The generalized distribution functions of the factorization approach could be used at temperatures 
up to about 10 20 K for such a system [fljif . 

All the above remarks led to recompute some bounds already obtained within the Tsallis et al. approach 
in this simpler framework, in order to get more reliable results and to check for consistency. This is what 
we briefly do below concerning the early universe test proposed in JO]. 

The test devised consist in compute a first order deviation in (q — 1) to the temperature of freezing out of 
the weak interctions in the early universe, Tf. This temperature is essential in the primordial nucleosynthesis 
scenario which, basically, is the competition between the rate of the expansion of the universe and that of 
the weak interactions which regulates the conversion of neutrons into protons and viceversa. When the 
expansion exceeds the rate of interactions they freeze out, and the final yields for the element production 
are roughly the ones we observe today. In particular, how a deviation in Tf cause a different prediction for 
the helium primordial yield, Y p , was made by Casas et. al. JET]. The result is, 
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Here, a radiation era relationship between time and temperature of the form (T cx t~? ) is assumed 
one sets 5T nuc — 0, because it is fixed by the binding energy of the deuteron. A = exp(— (t nuc — t 
stands for the fraction of neutrons which decayed into protons between tf and t nuc , with tf (t nuc ) the time 
of freeze out of the weak interactions (nucleosynthesis) and r the neutron mean lifetime. Considering now, 
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0.23 and \SY p \ = 0.01, which is the observational error for Y p 



values for the times and the mean life of neutron -which in fact, is not modified at order (q 
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A more detailed account of the processes that occurs in the early universe within this statistical framework 
is given elsewhere fl3| . 

To compute the {q — 1) corrections to Tf we recall the output of the factorization approach. The quantal 
distribution functions are given by, 



1 



[bosons] 



2 c x 



(1 — q) are 
2 (e x - l) 2 



(3) 



[fermions] 



l 



(1 - q) x 2 e x 
2 (e x + l) 2 



(4) 



There are two main corrections acting upon Tf. The first comes from the computation of the energy density 
of the universe. When the particles are higly relativistic, T ^> m, and non-degenerate T ^> fj,, we get 
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P fermions 
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gbj stands for the degeneracy factor of each one of the species involved. Using (||J^), we finally obtain 



Ptotal — Pbosons + Pf ermions 



30 



ff T 4 + 35.85 T 4 (q- 1), 



(7) 



where g = gt + 7/8<7/. At high enough temperatures, the energy density of the universe is essentially 
dominated by e~, e + , v and v and so gb = 2 and gf = 2 + 2 + 2 x 3 = 10. 



2 



The second correction comes from the computation of the weak interaction rates. We shall denote by 
Apn(T) the rate for the weak processes to convert protons into neutrons and by X np (T) the rate for the 
associated, reverse ones. Within the standard approximations applicable in the early universe regime 11 
it is possible to see that the weak interaction rate A is given by A ~ 4A v+n ^ p+e - . This last rate for the 
particular reaction quoted must be computed using, 



DC: 



K+n^p+e- = A j dp„p u p e E e (l- < n e >)n q (v) (8) 
Jo 

where A is a constant fixed by the experimental value of the neutron lifetime. Using (|3j^|), we get 

K + n^ P+ e- = K+nXte- + S 54 - 875 A (Q - 1), (9) 

as the biggest correction. Multiplying this result by 4, we obtain, 

^ = 1419.2 T 5 (1 - q). (10) 

Having in hands the main corrections that non-extensive quantal distribution functions provide, we move 
onwards to get the final bound upon q. To do so, we first compute the first order correction to Ty, which is 
defined as the temperature where the equality 



\ ( - ) . ^pt.. mi 



holds. The result is 

T- 



.,, -6.61(g-l). (12) 



Using, finally, Eq. (|2j) we get the following bound, 

\q- 1| < 4.01 x 10~ 3 . (13) 

The previous bound is, as one should expect, more restrictive than that obtained within the Tsallis et al. 
approach. [] This conclude the objective of this brief letter. 

Summing up, in this work we have revisited the recently proposed early Universe test of TT and computed 
the related bound on q using the generalized distribution functions within the factorization approximation. 
Our main result is written in Eq. (|i~3|): 1 to 100 seconds after the Big Bang, q must be that close to 1 
in order to be able to reproduce observational results on helium abundance. As it is expected, the bound 
is found to be consistent with those of Tsallis et al. approach, a result which was also obtained in Rcfs. 
|js|,fio| . Although all these efforts clarify the fact that the factorization approach results can be used with 
confidence for physical systems, our belief is that it will be used much more in the near future especially for 
the applications with q values far from unity, where Tsallis et al. approach is unapplicable. Therefore new 
attempts on this line would be highly welcomed.^] 
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1 Recall the Errata in jTT]. The final bound turns out to be |1 — q\ < 3.4 x 10" 3 . 

2 After submission of this work we became aware of the work by Rajagopal et al. [Phys. Rev. Lett. 80, 3907 (1998)] 
where exact results for quantal distribution functions are given. We have explored how these results match with 
Tsallis et al.'s and factorization approaches and we shall report on it in a future comunication. 
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